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Abstract 

Background: Sexual size dimorphism (SSD) is widespread and variable among animals. Sexual selection, fecundity 
selection and ecological divergence between males and females are the major evolutionary forces of SSD. However, 
the influences of mating system and habitat types on SSD have received little attention. Here, using phylogenetic 
comparative methods, we at first examine the hypotheses to that mating system (intensity of sexual selection) and 
habitat types affect significantly variation in SSD in anurans (39 species and 18 genera). 

Results: Our data set encompass 39 species with female-biased SSD. We provide evidence that the effects of mat- 
ing system and habitat types on SSD were non-significant across species, also when the analyses were phylogenet- 
ically corrected. 

Conclusions: Contrast to the hypotheses, our findings suggest that mating system and habitat types do not play 
an important role in shaping macro-evolutionary patterns of SSD in anurans. Mating system and habitat types can- 
not explain the variation in SSD when correcting for phylogenetic effects. 
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Introduction 

Sexual size dimorphism (SSD) is widespread and variable 
among animals [1]. In some groups (e.g. birds, lizards 
and most mammals), males are bigger than females, 
whereas in other groups (fishes and anurans) females are 
bigger than males [2], Some key hypotheses are provided 
for explaining the evolution and maintenance of SSD, al- 
though their explanatory power remains controversial 
[3-5]. It is now widely agreed that sexual selection in 
favor of large males to improve intra-sexual combat suc- 
cess and fecundity selection for large females to increase 
reproductive output are the major evolutionary forces of 
SSD in many organisms [1]. Beyond the two hypotheses, 
ecological divergence between the sexes due to intraspe- 
cific competition has been proposed to explain evolution 
of SSD [6]. 

SSD is often used as an indicator of the intensity of 
sexual selection in animals [7]. Mating system is associ- 
ated with the intensity of sexual selection [1], The idea is 
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that promiscuous species where females mate with more 
than one male have more intensive competition than 
monogamous species where females mate with only one 
male. Sexual selection hypothesis predicts that the inten- 
sity of selection promotes variation of SSD among spe- 
cies through intra-sexual competition or inter-sexual 
mate choice favoring large size in one sex [1,7]. Conse- 
quently, mating system can mediate evolution of SSD. 
Moreover, habitats types have been proposed to affect 
SSD due to energy constraints and predators. For ex- 
ample, small male size in aquatic habitats resulting from 
selection to reduce energy expenditure in mate searching 
shows a female-biased SSD, whereas large male size in 
terrestrial habitats results from predation pressure, mate 
searching needs, or desiccation avoidance, showing a 
mixed SSD [6,8,9]. For frogs it appears that aquatic and 
arboreal species have usually smaller variation in SSD 
than terrestrial species [10]. 

Anurans inhabit a wide range of habitats (i.e. wetlands, 
grasslands, steams, trees and ponds) and exhibit a re- 
markable diversity of mating system (i.e. social polyandry 
and social monogamy) that is unique among vertebrates 
[10,11]. If the selection hypotheses are valid, we expect 
mating system and habitat types having effects on SSD 
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in anurans. Comparative studies on anuran SSD have 
been shown that differences in the age structure between 
the sexes in breeding populations can explain variation 
in SSD [12-15]. Moreover, several attempts have been 
made to explain SSD as a consequence of sexual selec- 
tion, fecundity selection and life-history traits in anurans 
[16,17]. However, the influences of mating system and 
habitat types on SSD in anurans have received little at- 
tention. Here, we at first examine the hypotheses to that 
mating system and habitat types significantly influence 
on variation of SSD. 

Results 

Thirty-nine species were characterised by female-biased 
SSD. The GLM revealed that the mean SVL significantly 
differed between the sexes (F h 77 = 3.059, P = 0.043) and 
terrestrial, arboreal, semiaquatic and aquatic habitats 
(^3, 77 = 2.756, P = 0.050), but did not differ among sequen- 
tial polyandry, simultaneous polyandry and monoandry 
(^2, 77 = 0.641, P = 0.530). However, there were also non- 
significant mating systems*sex interaction (F 2) 77 = 0.036, 
P= 0.964) and habitats *sex interaction (F 3f 77 = 0.125, 
P = 0.945), revealing that the degree - but not direction (fe- 
males always the larger sex) - of SSD did not differ among 
mating system and habitats (Figure 1 and 2). 

We tested for the effects of mating system and habitats 
on SSD contrasts using generalized least squares. We 
found the mean SVL contrasts significantly differed be- 
tween the sexes (F h 75 = 3.182, P = 0.041) and habitats 
(F 3) 75 = 4.342, P = 0.020), but did not differ among mat- 
ing system (F 2 , 75 = 0.038, P = 0.963). However, the non- 
significant interactions between mating system and sex 
(F2, 75 = 0.017, P = 0.983) and between habitats and sex 
(F 3} 75 = 0.042, P = 0.997) revealed the evidence that in- 
tensity of sexual selection and habitats were not associ- 
ated with variation of SSD contrasts. 
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Figure 1 Sexual dimorphism in body size [mean±SD log 10 
(female size)-log 10 (male size)] in relation to mating system 

among 39 anuran species. 
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Figure 2 Sexual dimorphism in body size [mean ± SD log 10 
(female size)-log 10 (male size)] in relation to habitat types 
among 39 anuran species. 1 = Terrestrial - mostly occur on ground, 
forage in ground in various conditions, 2 = Semiaquatic - not entirely 
aquatic, usually living or growing in or near water, 3 = Arboreal - 
mostly occur on trees, forage in trees and rarely come down to the 
ground, 4 = Aquatic - mostly occur on water, forage in water 
in stream. 
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Discussion 

Ninety percent anuran species is characterised by 
female-biased SSD, and 10% by male-biased SSD [2]. In 
our study, female-biased SSD is the predominant pat- 
tern. Sexual selection supports the idea that males en- 
gaging in physical combat with one another may select 
for large males when such males are more successful in 
gaining access to mates [1]. However, male-male com- 
petition prevails in some anurans where a female-biased 
SSD is common [16]. This means that selection on large 
females should be stronger than selection on large males 
due to the expression of a given SSD depending on rela- 
tive dominance of competing selective forces [17]. Such 
a source of selection in favor of large females is consid- 
ered to be the size-dependent fecundity advantage, 
which can drive the evolution of female body size and 
consequently leads to the evolution of female-biased 
SSD [16]. Moreover, females begin breeding later, live 
longer but grow more slowly than males, resulting in 
female-larger patterns of SSD in anuran lineages. In par- 
ticular, SSD is increasingly biased towards females across 
species when the duration of growth in females is longer 
than males [18]. 

Mating system is often used as an indicator of the in- 
tensity of sexual selection [1]. Our results show that SSD 
is not affected by mating system, suggesting that anurans 
with intense sexual selection do not exhibit a range of 
dimorphisms. This pattern is contradictory to the well- 
established concept that male-male competition (sexual 
selection) drives to increase male body size and results 
in male-biased SSD [1]. Several potential causes can ex- 
plain why intensity of sexual selection does not promote 
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evolution of SSD in anurans [1,11,15,19]. Firstly, mating 
success is positively correlated with male body size in 
some species, but not in others. Secondly, an advantage 
for small males relates to post-copulatory sexual selec- 
tion through sperm competition. The idea is that small 
males have low chances of obtaining females through fe- 
male choice or male-male competition, so they should 
instead disproportionately invest in mate search and 
sperm competition, where they are assumed to have 
relatively better chances. Thirdly, the occurrence of al- 
ternative mating tactics in anurans reduces the selection 
pressure on male body size. Consequently, small males 
may gain fitness by using alternative mating tactics. Fi- 
nally, the extreme diversity of life-history traits and their 
plasticity may mask the potential contribution of male- 
male competition to body size evolution. 

Mating system mediates SSD in birds, such that under 
polyandry, sexual selection on females results in a SSD pat- 
tern opposite to Renschs rule while selection on males re- 
sults in a SSD pattern consistent with Renschs rule [20]. 
Evidence in support of this idea has been obtained from 
Renschs rule is driven by a correlated evolutionary re- 
sponse in one sex to stronger size selection in the other 
sex. However, contrast to the hypothesis, we find that mat- 
ing system in anurans cannot mediate SSD, such that sex- 
ual selection on males results in a SSD pattern inconsistent 
with Renschs rule and its reverse [13]. The possible reason 
is that fecundity selection for large females balances out 
sexual selection on large males. 

Anurans are a diverse group of vertebrates renowned for 
variable life-history traits, which include mainly terrestrial, 
arboreal, semiaquatic and aquatic habitats [10]. For ex- 
ample, terrestrial species exhibit male combat, and males 
are as large or larger than females. For aquatic and arboreal 
species that exhibit female choice due to male mating calls, 
males are smaller than females. Consequently, terrestrial 
species should show larger variation in SSD than aquatic 
and arboreal species [10,21]. However, our results suggest 
that in anurans the regimes of natural selection imposed by 
habitat types alone may not have exerted a significant im- 
pact on body size in either of the sexes. This finding sup- 
ports the assumption that males and females are 
ecologically or phenotypically equivalent which may pro- 
vide an incomplete or even mistaken picture of the process 
of body size diversification [22]. In conclusion, variation of 
SSD in anurans cannot be explained by mating system and 
habitat types. It is the result of a variety of selective forces, 
including sexual selection, fecundity selection, life-history 
and ecological factors. 

In this study, methodological aspects give reason to view 
our results with caution. At first, the phylogenetic tree we 
present appears to be a simple dendrogram showing just 
branching pattern, i.e. its branch spans do not represent 
time or relative amount of character change. PGLM 



methods make explicit use of information contained in 
branch lengths. However, we cannot obtain a tree with 
branch lengths, so we use GLM to test variation of SSD 
contrasts which may outcome potential problems and 
biases which may occur as a consequence of ignoring this 
information. At second, correlations across species should 
be regarded cautiously based on the fact that species' data 
points cannot be assumed to be statistically independent 
[23,24]. However, comparisons across species still result in 
meaningful analyses unless they need a cluster of points 
that share an immediate common ancestor Harvey and 
Pagel [25]. 

Materials and methods 

We obtained sex-specific demographic and morpho- 
logical data on mean size across 39 species and 18 gen- 
era from the literatures (Additional file 1: Table SI). We 
calculated the mean values for the population as alge- 
braic means for each year, weighted by sample size. 
Mean values for species were obtained as algebraic 
means of population values regardless of the sample size 
in cases where data were available for different popula- 
tion [12]. Following the method proposed by Roberts 
and Byrne [26], we used mating system as an imperfect 
surrogate of the intensity of sexual selection on a three- 
point scale: 1 = sequential polyandry where two or more 
males simultaneously releasing sperm or sequentially re- 
leasing sperm in a time frame that allows for the occur- 
rence of sperm competition; 2 = simultaneous polyandry 
where a females mates with two or more males over the 
course of a breeding season by depositing part of a sin- 
gle clutch with each male or, multiple clutching; 3 = 
monoandry where a females mates with one male over 
the course of a breeding season by depositing part of a 
single clutch. Habitat types were classified on a four- 
point scale: 1 = Terrestrial - mostly occur on ground, 
forage in ground in various conditions, 2 = Semiaquatic - 
not entirely aquatic, usually living or growing in or near 
water, 3 = Arboreal - mostly occur on trees, forage in 
trees and rarely come down to the ground, 4 = Aquatic 
- mostly occur on water, forage in water in stream. Fol- 
lowing the methods by Lovich and Gibbons [27], we cal- 
culated SSD as (log 10 (female mean size)/log 10 (male 
mean size)) - 1, arbitrarily set positive when females are 
larger and negative when males are larger. All animals 
used in this study were treated humanely and ethically 
following all applicable institutional Animal Care guide- 
lines in China. 

Comparative analyses of interspecific data may require 
phylogenetic control as closely related species share parts of 
their evolutionary history. Therefore, they cannot be consid- 
ered independent data points for statistical analyses [24]. 
Phylogenetic analyses were based on generalized least 
squares, which is a powerful and comprehensive approach 
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to the analysis of comparative data [28]. Generalized least 
squares is a modification of generalized linear models 
(GLMs) in which the phylogeny is used to specify the ex- 
pected variance and covariance between species under an 
assumed. For our comparative analysis, we used an estab- 
lished phylogeny [29,30] (Additional file 2: Figure SI). We 
calculated mean size in both sexes for ancestral nodes as 
the algebraic mean of the two closest lower nodes [23]. De- 
tails of the general procedure for estimating the character 
values in the ancestors are presented in Felsenstein [24]. 
With 39 species at the tips of this reconstructed tree, 38 
(39-1) body size of contrasts within each sex could be 
computed for pairs of nodes sharing an immediate com- 
mon ancestor, and then re-scaled and analysed as suggested 
by Garland et al. [31]. Correct standardization and homo- 
geneity of variance of standardized contrasts were con- 
firmed using the method proposed by Purvis and Rambaut 
[32]. To conduct the effect of mating system and habitat 
types on variation in SSD, we conducted conventional non- 
phylogenetic GLMs separately on log (body size) as a 
dependent variable, habitats, mating system, sex and their 
interactions as fixed factors. We then used GLMs separ- 
ately on log (body size) contrasts as a dependent variable, 
and mating system, habitats, sex and their interactions as 
fixed factors to test the effects of mating system and habitat 
types on SSD contrasts. All tests were conducted by using 
Type III sums of squares. 

Additional files 



Additional file 1: Table SI. Species, sample size (females/males), mean 
body size in males and females of 39 species and references of published 
papers. 

Additional file 2: Figure SI. The phylogenetic tree of the 39 anuran 
species used in the comparative analysis following Frost et al. (2006) and 
Pyron and Wiens (2011). 
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